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HIGHLIGHTS 


►  A  life  model  is  developed  by  including  the  Mn  dissolution  for  a  lithium  ion  cell  with  spinel  cathode. 

►  The  radius  change  due  the  Mn  dissolution  is  considered  in  the  life  model. 

►  The  temperature  effect  on  the  Mn  dissolution  is  included  in  the  life  model. 

►  The  resistance  of  the  inactive  film  generated  due  to  the  Mn  dissolution  is  included  in  the  life  model. 
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A  life  model  is  developed  for  a  lithium  ion  cell  with  a  spinel-based  cathode.  It  is  assumed  in  the  proposed 
model  that  the  Mn(III)  disproportionation  reaction  causes  the  degradation  of  the  spinel  cathode.  The 
Mn(III)  disproportionation  reaction  leads  to  the  Mn(II)  dissolving  into  the  electrolyte  and  the  formation 
of  an  inactive  material  layer  which  causes  a  resistance  increase  in  the  cathode.  The  proposed  model  is 
used  to  investigate  the  effects  of  ambient  temperature  and  voltage  range  of  cycling  on  the  loss  of  the  cell 
capacity  and  the  changes  in  the  volume  fraction  of  the  cathode  active  material,  the  radius  of  the  cathode 
particle  and  the  resistance  of  the  cathode. 

©  2012  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

Spinel  LiMn2C>4  is  of  great  interest  as  a  cathode  material  for 
rechargeable  lithium-ion  batteries  because  of  its  low  cost,  high 
potential  and  safety.  Especially,  spinel  LiMn2C>4  maybe  play  great 
role  in  the  near  future  for  electric  vehicle  applications.  However, 
this  material  exhibits  serious  capacity  fading  during  cycling  or 
storage  at  elevated  temperature  which  makes  it  less  competitive 
with  other  cathode  materials. 

Recent  studies  reveal  that  the  capacity  fading  in  spinel  LiMn204 
cathode  is  attributed  to  manganese  dissolution  [1,2],  There  are  at 
least  two  possible  reasons  for  Mn  dissolution:  (1 )  Mn(II)  dissolution 
into  the  electrolyte  due  to  acid  attack  and  a  disproportion  reaction 
of  Mn(III)  at  the  particle  surface  [3];  (2)  the  instability  of  the  two- 
phase  structure  in  the  charged  state  which  leads  to  the  loss  of 
MnO  and  the  dissolution  of  Mn(II)  to  form  a  more  stable  single¬ 
phase  structure  [4],  In  the  cycle  performance  studies  on  the 
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spinel  LiMn204  cathode,  Jang  et  al.  [5]  and  Xia  et  al.  [4]  found  that 
manganese  dissolution  was  notably  high  at  the  charged  state 
(delithiated  state)  of  the  spinel  LiMn2C>4  cathode.  Wang  et  al.  [6] 
used  rotating  ring-disk  collection  experiments  and  found  that  the 
concentration  of  dissolved  Mn(II)  increases  with  cycle  number  and 
elevated  temperature.  In  any  cycle,  the  ring  currents  exhibited 
maxima  corresponding  to  the  end-of-charge  (EOC)  and  end-of- 
discharge  (EOD),  while  the  largest  peak  occurred  at  the  EOC  state. 
The  Mn  dissolution  from  the  spinel  LiMn204  cathode  occurs  during 
cycling,  especially  in  the  charged  state  (i.e.,  the  cathode  potential  is 
greater  than  4.1  V)  and  in  the  discharged  state  (i.e.,  the  cathode 
potential  is  less  than  3.1  V). 

The  storage  performance  of  a  spinel  LiMn204  cathode  was  also 
studied.  Aoshima  et  al.  [7]  investigated  the  capacity  loss  of  the 
spinel  cathode  at  various  states  of  charge  (SOC)  storage  at  70  °C.  The 
most  significant  capacity  loss  was  observed  in  the  case  of  the  dis¬ 
charged  cathode.  Inoue  and  Sano  [8]  also  found  similar  results. 
Their  storage  experiments  at  80  °C  with  different  additives  revealed 
that  water  can  accelerate  the  capacity  fading  of  the  spinel  cathode. 

Electrochemical  impedance  spectroscopy  (EIS)  studies  revealed 
that  the  rate  capacity  loss  due  to  the  resistance  increase  is  another 
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important  contribution  due  to  spinel  degradation  [4,5,9—12].  The 
resistance  increase  in  the  spinel  cathode  resulted  from:  the  SE1  film 
building  up  due  to  the  oxidation  of  electrolyte  and  the  loss  contact 
with  the  conducting  additive  due  to  the  formation  of  inactive 
material  in  the  Mn(III)  disproportionation  reaction. 

A  mathematical  model  was  developed  by  Park  et  al.  [13]  to 
simulate  the  degradation  of  spinel  LiMn2C>4  cathode  based  on  the 
mechanism  of  Mn(III)  disproportionation  reaction  by  Lu  and  Lin 
[14]  Although  the  volume  change  in  the  active  spinel  cathode  was 
taken  into  account,  the  change  in  the  porosity  of  the  spinel  cathode 
due  to  the  Mn  dissolution  was  not  included  in  their  model. 

A  comprehensive  pseudo  two-dimensional  (P2D)  electro¬ 
chemical  thermal  model  for  the  degradation  of  the  spinel  LiMn2C>4 
cathode  due  to  Mn(IIl)  disproportionation  reaction  is  developed  in 
this  paper.  The  proposed  model  includes  the  volume  changes  in 
both  the  active  material  and  the  inactive  material  in  the  cathode, 
the  radius  change  in  the  cathode  solid  particles,  the  porosity  change 
in  the  cathode,  and  film  resistance  increase  in  the  cathode  due  to 
Mn(III)  disproportionation  reaction.  The  proposed  model  was 
applied  to  simulate  the  cycling  performance  at  different  tempera¬ 
tures  and  voltage  ranges  for  a  lithium-ion  cell  which  consists  of 
a  spinel  LiMn204  cathode,  Celgard’s  25  pm  (in  thickness)  separator, 
a  carbon-based  anode  and  1  M  LiFP6  in  PC/EC/DMC  (1:1:1  volume 
ratio)  electrolyte. 

2.  Model  development 

2.1.  Modeling  ofMn  dissolution 

ft  is  assumed  that  the  Mn  dissolution  occurs  according  to  the 
following  disproportionation  reaction  of  Mn(III)  [5,13,14]: 

2Mn(III)(s)  — >Mn(II)(1)  +  Mn(IV)(s)  (1) 

where  Mn(III)(s)  is  the  active  Mn  with  the  valence  of  3  in  the  spinel 
cathode,  Mn(Il)(i)  is  the  product  in  the  disproportionation  reaction 
with  the  valence  of  2  and  is  soluble  in  the  electrolyte,  and  Mn(IV)(S) 
is  another  product  with  the  valence  of  4  and  remains  in  the  solid 
phase  but  is  inactive  for  the  Li+  insertion/deintercalation  reaction. 
The  disproportionation  reaction  of  Mn(III)  is  not  an  electrochemical 
reaction  and  occurs  all  the  time  (except  for  a  fully  charged  cell 
where  all  Mn  is  Mn(lV))  as  long  as  the  cell  is  assembled.  The 
reaction  rate  of  the  disproportionation  reaction  depends  on  the 
temperature  only.  It  was  proposed  that  the  conversion  of  reaction 
(1)  is  a  function  of  time  and  is  given  by  [14]: 

1  -  (1  -Xa)1/3=  kt  (2) 

where  k  is  the  reaction  rate  constant  (1/s),  andXa  is  the  conversion 
of  the  dissolution  reaction  shown  in  Eq.  (1)  and  is  defined  by: 


The  cathode  active  material,  LiMn2C>4,  can  be  rewritten  as  LiM- 
n(III)Mn(IV)C>4,  where  Mn(III)  is  the  active  part  which  reacts  with  Li 
in  the  electrochemical  reaction.  Rewrite  Eq.  (1)  in  terms  of  the 
active  material  and  the  products  with  the  same  number  of  Mn 
atoms  in  the  dissolution  reaction  as: 

4A— >3S  +  L  (5) 

where  A  denotes  the  active  material  (LiMn2C>4)  which  is  consumed 
due  to  the  dissolution  reaction,  S  is  the  solid  phase  product  with  two 
Mn  atoms  in  4  valence  (Mn(IV)2)  and  L  is  the  liquid  phase  product 
with  two  Mn  atoms  in  2  valence  (Mn(II)2).  Because  A,  S  and  L  consist 
of  two  Mn  atoms,  it  is  assumed  that  A,  S  and  L  have  approximately 
the  same  molar  volumes.  We  denote  the  remaining  active  material 
(LiMn2C>4)  as  R.  Thus,  the  sum  of  volumes  of  A  and  R  is  a  constant 
which  is  the  initial  volume  of  the  active  material:  V0,  that  is: 

V0  =  VR  +  VA  (6) 

where  VR  and  VA  denote  the  volumes  of  R  and  A,  respectively. 
According  to  the  stoichiometry  in  Eq.  (5)  and  along  with  the 
assumption  that  the  molar  volumes  of  A,  S  and  L  are  the  same,  we 
have  the  following  relationships: 

Vi  =  lvA  (7) 


Eq.  (3)  can  be  rewritten  as: 


(2  moles  L)  x  2  _  4  mole  L 
3  ~  1  mole  R  —  1  mole  R 

Since  we  assumed  that  the  molar  volumes  of  L  and  R  are  the 
same,  Eq.  (9)  is  equivalent  to  the  following  equation: 


Xa 


4Vl 

Vr 


Substitution  of  Eq.  (7) 


Eq.  (10)  yields: 


(10) 


VA 

Vr 


(11) 


Eq.  (11)  indicates  that  the  conversion  of  Mn  dissolution  reaction 
equals  the  volume  of  the  dissolved  active  material  divided  by  the 
volume  of  the  remaining  active  material. 

The  initial  volume  fraction  of  the  active  material  in  the  spinel 
cathode  is  defined  by: 


(Number  of  moles  of  the  dissolved  Mn(II)  in  the  electrolyte)  x  2 
a  —  Number  of  moles  of  the  Mn(lll)  in  the  spinel  cathode 


(3) 


The  reaction  rate  constant  is  defined  by: 


(4) 


(12) 


where  k0  is  the  pre-exponential  coefficient,  s  and  Ear  is  the 
activation  energy  related  to  k. 


where  Vcath  is  the  volume  of  the  cathode.  The  volume  fraction  of  the 
remaining  active  material  is  given  by: 
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= 

>cath 

By  substituting  Eqs.  (7),  (8),  (11 )  and  (12)  into  Eq.  (13), 


(13) 


The  volume  fraction  of  the  inactive  material,  qv,p, 
phase  is  defined  by: 


(14) 
n  the  solid 


**  =  V 

The  porosity  of  the  cathode, 


where  Rp  is  the  initial  radius  of  the  particle  and  V0  =  4/3tcR^. 
Based  on  Eq.  (20),  the  radius  of  the  solid  particle,  Rs,p(t),  as  a  func¬ 
tion  of  time  is  given  by: 


Rs 


s,P  (t)  =  ( 


1  +  0.75Xa\ 
1  +  Xa ) 


1/3 

RP 


(22) 


Since  the  inactive  solid  phase  grows  as  the  dissolution  reaction 
occurs,  the  resistance  of  this  inactive  material  layer  increases. 
Consequently,  we  introduce  a  variable,  Rfiim,  to  denote  the  film 
resistance  in  the  cathode  particle.  We  also  assume  that  this  film 
resistance  depends  on  Xa  as  follows: 


Rfiim  =  Rfllm  +  RSEI 


Rs,p(f)  —  Ra,p(f) 

Rp 


(23) 


£f,P 


(16) 


where  tfP  is  the  volume  fraction  of  the  filler  in  the  cathode.  By 
substituting  Eqs.  (14)  and  (15)  into  Eq.  (16),  we  have: 


£P 


o  1  +  0.75Xa 
1  +Xa 


(17) 


The  radius  of  the  active  material  in  the  cathode  solid  particle 
and  the  radius  of  the  solid  particle  are  illustrated  in  Fig.  1. 

The  volume  of  the  remaining  active  material  as  a  function  of  the 
initial  volume  of  the  solid  particle  is  obtained  by  substituting  Eq. 
(11)  into  Eq.  (6): 


Vr 


Vp 

(1+Xa) 


(18) 


The  volume  of  the  solid  phase  as  a  function  of  time  is  given  by: 


V(t)  =VS  +  VR  (19) 

Substitution  of  Eqs.  (7),  (11)  and  (18)  into  Eq.  (19)  yields: 


V(t)  = 


1  +  0.75Xa 
1  T  Xa 


Vo 


(20) 


According  to  Eq.  (18),  the  relationship  between  the  radius  of  the 
active  particle  and  the  initial  radius  of  the  particle  is  given  by: 


/  1  \V3 
Ra’P(t)  =  (iTXa) 


(21) 


where  Rfj|m  is  the  initial  film  resistance  and  Rsei  is  the  resistance  of 
the  inactive  material. 

2.2.  Electrochemical  model 

A  schematic  of  a  typical  Li-ion  cell  consisting  of  three  regions 
namely  a  positive  electrode  (spinel  LiM^CU),  a  separator  (poly¬ 
propylene),  and  a  negative  electrode  (carbon-based)  is  shown  in 
Fig.  2.  A  physics-based  pseudo  two-dimensional  (P2D)  electro¬ 
chemical  model  which  developed  by  Newman  et  al.  [15]  was 
applied  to  simulate  a  lithium  ion  cell  with  a  spinel  cathode  and 
a  carbon-based  anode  as  shown  in  Fig.  2.  The  P2D  electrochemical 
model  consists  of  the  material  balance  in  the  solid  particles  in  the 
micro  scale  (r-direction  in  Fig.  2),  the  material  balance  in  the  binary 
electrolyte  which  fills  the  pores  in  the  porous  electrodes  and  the 
separator,  and  the  charge  balances  in  the  solid  matrix  in  the  elec¬ 
trodes  and  in  the  electrolyte.  The  governing  equations  of  the  P2D 
model  are  summarized  in  Appendix  A  besides  of  the  material 
balance  of  Li+  in  the  solid  particles  in  the  cathode.  The  open  circuit 
potentials  of  the  spinel  cathode  and  the  carbon-based  anode  and 
other  model  expressions  are  summarized  in  Appendix  B. 

The  simulation  of  lithium  ions  diffusion  in  the  cathode  material 
with  multi-phases  in  a  lithium  ion  cell  has  been  widely  investigated 
by  using  the  sharp  interface  models  (i.e.,  the  shrinking  core  model 
[16]  and  the  moving  boundary  model  [17])  and  the  phase  field 
model  [18]  which  is  a  non-Fickian  behavior  and  interfaces  tracking 
free  model.  We  assumed  that  the  Mn  dissolution  reaction  (Reaction 
(1))  only  occurs  on  the  interface  between  the  active  spinel  and  the 


Fig.  1.  Radius  change  in  a  solid  particle  due  to  Mn  dissolution. 


Fig.  2.  Schematic  of  a  Li-ion  cell. 
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electrolyte.  This  produces  an  inactive  solid  porous  layer  (shown  in 
Fig.  1 )  on  the  surface  of  the  particle.  Consequently,  the  electrolyte 
can  penetrate  the  inactive  porous  layer  and  reach  the  surface  of  the 
active  spinel  where  the  Li  insertion/extraction  reaction  occurs. 
Based  on  these  assumptions  a  moving  boundary  model  is  applied  to 
simulate  the  diffusion  of  Li  ions  in  the  active  spinel.  The  material 
balance  of  Li+  in  the  solid  particles  in  the  cathode  with  a  moving 
reaction  boundary  which  is  shown  in  Fig.  1  is  given  by: 


with  the  boundary  conditions: 


Qrev.i  —  ftfJjT; 


a 


(32) 


where  [at//9T]f  is  the  open  circuit  potential  change  rate  due  to  the 
change  in  temperature  in  electrode  i  evaluated  at  the  reference 
temperature  Tref.  The  expressions  of  [011/01];  as  functions  of  state  of 
charge  for  the  electrodes  are  presented  in  Appendix  B. 


1  OCj  d<p2.i  . 


-Ds,p%a|  =  0  (25) 

(,r  lr=0 

-0s,p%1  =jp  (26) 

O'  lr=Ra,p(t) 

and  the  initial  condition: 

cs,p(t  =  0)  =  cSiP,o  (27) 

where  Ds  p  is  the  diffusion  coefficient  of  Li+  in  the  solid  particles  in 
the  cathode,  jp  is  the  pore  wall  flux  in  the  cathode,  and  cs,p,o  is  the 
initial  concentration  of  Li+  in  the  particles. 

Since  the  radius  of  the  active  part  of  the  cathode,  Ra,p,  is  varying 
with  time,  it  is  more  convenient  to  convert  Eqs.  (24)— (27)  into  the 
dimensionless  spatial  coordinate,  y,  which  is  defined  by: 

y  =  ^  (28) 

«a,p 

Eqs.  (24)— (26)  in  the  dimensionless  spatial  coordinate,  y,  with 
fixed  boundaries  are  given  by: 


The  initial  condition  for  the  temperature  is  given  by  the  ambient 
temperature: 


where  Tair  is  the  ambient  temperature.  The  boundary  conditions  are 
determined  by  Newton’s  cooling  law: 


-v£|  =  hcoef(T  -  Tair)  (37) 

dXlx=Lp+Ls+Ln 

where  hcoef  is  the  heat  transfer  coefficient. 

Table  1  shows  the  values  of  the  model  parameters. 


0CSiP 

9t 

~Ds,  P 


1  8 

Rfp Fey 


0Cs,p| 

’^ly=0 


=  0 


_Ds’P“air  .  -  Ravip 


(29) 


2.3.  Thermal  model 

The  energy  balance  for  the  cell  is  given  by  [19]: 

PfCp.ff  =  ^  (/igx)  +  Qrxn’’  +  +  Qohnu  (30) 

where  pi  is  the  density  of  cell  component  i  (i  =  p  for  the  cathode, 
i  =  s  for  the  separator,  and  i  =  n  for  the  anode),  CPij  is  the  specific 
heat  of  the  component  i,  A;  is  the  thermal  conductivity  of  compo¬ 
nent  i,  Qrxn,j  is  the  heat  generation  rate  due  to  electrochemical 
reaction  in  electrode  z,  Qrev.i  is  the  heat  generation  rate  due  to 
reversible  reaction  in  electrode  i,  and  Qohm,i  is  the  heat  generation 
rate  due  to  Ohmic  loss  in  component  i.  The  three  types  of  heat 
generation  rates  are  defined  as  follows: 


Table  1 

Model  parameters. 


Parameter  Value  Parameter  Value 


h. 

u 

Ln 

V 

Rna 


4” 

4" 

R 


bruggpb 

bruggsb 

bruggnb 


135e— 6  m  <rpa 

25e— 6  m  <r„a 

100e-6  m  Ds,p,refa 

Se-6  m  Ds.„,refa 

12.5e— 6  m  kprJ’ 

0.444  knrefb 

0.41  koc 

0.357  Earc 

0.252  Pad,  Pb 

0.172  £ad,nb 

1  %  f/.P  Lik.p1’ 


8.314  J  mar1  IT1  /sd 

964,87  C  equiv-1  V* 

298.15  K  ppd 

328.15 1<  psd 

0.363  />„“ 

1.5  Cp,pd 

1.5  Cp/ 

1.5  Cp’nd 

35  A  m  2  ai  v 

228,60  mol  m  3  atp 

263,90  mol  m  3  ain 

0.30  acn 


3.41e5  s1 
72.48e3  J  mor1 
2.9e3  J  mor1 
3.5e3  J  mor1 
5.8e3  J  mol-1 
2.0e3  J  mol-1 
1.58  W  m-1  K1 
0.344  W  m1  K1 
1.04  Wm1  K-1 
2328.5  kg  m  3 
1008.98  kg  irr3 
1347.33  kg  rrr3 
1269.21  J  kg1  !<-' 
1978.16  J  kg1  K_1 
1437.4  J  kg1  K-1 
0.5 
0.5 
0.5 
0.5 


0.75  Rglmb  1.0e-3  0  m2 

2000  mol  rrr3  hCOef  2.0  W  m  2  1C1 


a  Obtained  from  Ref.  [20], 
b  Assumed. 

c  Obtained  from  Ref.  [14], 
d  Obtained  from  COMSOL  model  library. 


Qrxn.i  =  FaJi  (01  i  -  <l>2.,  ~  Ui ) ,  i  =  P,  n 

where  Ui  is  the  open  circuit  potential  of  electrode  i. 


(31) 
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3.  Results  and  discussion 

The  IDA  package,  which  can  be  used  to  solve  initial  value 
problems  for  differential/algebraic  equation  systems,  in  SUNDIALS 
(Suite  of  Nonlinear  and  Differential/Algebraic  Equation  Solvers) 
[21  ]  was  applied  to  solve  the  P2D  thermal  model  associated  with 
the  proposed  Mn  dissolution  model.  The  finite  volume  method 
(FVM)  was  applied  to  discretize  the  spatial  coordinates.  By  using 
FVM,  the  cathode  and  anode  were  portioned  into  50  finite  volumes 
and  the  separator  and  the  particles  in  r-direction  (shown  in  Fig.  2) 
are  divided  into  25  finite  volumes.  The  SI  unit  system  was  applied  in 
the  system.  The  relative  and  the  absolute  tolerances  were  selected 
as  1.0e-9  for  all  of  the  dependent  variables  in  the  system.  The  in- 
house  C  code  was  validated  by  comparing  the  results  obtained 
from  the  in-house  C  code  to  those  obtained  from  COMSOL  [22]  in 
which  the  same  equations  were  solved  and  the  same  tolerances 
were  used. 

Four  case  studies  were  carried  out  using  our  in-house  C  code. 
The  cell  mentioned  in  the  previous  section  was  discharged  and 
charged  fifty  times  at  constant  current  rate  of  2C  (1C  =  17.5  A  m-2) 
with  two  different  ambient  temperatures:  25  °C  and  55  °C  and  two 
different  voltage  ranges:  [3.2,  4.0]  V  and  [3.5,  4.3]  V. 

Fig.  3  shows  the  normalized  capacity  obtained  during  the  2C 
discharge  vs.  the  cycle  number  for  the  four  cases.  The  normalized 
capacity  in  the  ith  cycle,  NCap,-,  is  obtained  as  follows: 

NCap,,ga  (38) 

where  Neap,-  is  the  normalized  discharge  capacity  in  cycle  i,  Cap,-  is 
the  discharge  capacity  in  cycle  i  and  Capi  is  the  discharge  capacity  in 
cycle  1.  Fig.  3  indicates  that  the  cell  capacity  loss  is  more  significant 
due  to  the  spinel  cathode  degradation  when  the  cell  is  cycled  at 
a  higher  ambient  temperature  within  the  same  voltage  range.  It  is 
also  found  from  Fig.  3  that  at  the  same  ambient  temperature,  the  cell 
cycled  in  a  higher  voltage  range  (i.e.,  3.5— 4.3  V)  loses  more  capacity. 
As  mentioned  in  the  previous  section,  we  assumed  that  the  Mn 
dissolution  reaction  does  not  depend  on  the  state  of  charge  of  the 
spinel  cathode.  The  difference  in  the  voltage  range  only  affects  the 
discharge  time  and  the  cell  temperature  which  cause  the  difference 
in  the  degradation  rate  of  the  spinel  cathode.  The  discharge  time  for 
a  cell  discharged  at  2C  in  the  voltage  range  of  3.5-4.3  V  is  longer 
than  that  fora  cell  discharged  at  2C  in  the  voltage  range  of  3. 2-4.0  V. 
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Fig.  4.  Averaged  volume  fraction  of  the  active  cathode  at  the  end  of  discharge  vs.  the 


Also  the  cell  temperature  is  a  slightly  higher  in  the  cell  cycled  within 
3.5— 4.3  V  than  in  the  cell  cycled  within  3.2-4.0  V.  The  capacity  loss 
shown  in  Fig.  3  is  the  combination  of  the  contributions  of  the  loss  of 
the  active  spinel  cathode  due  to  Mn  dissolution  reaction  and  the 
resistance  increase  on  the  spinel  cathode. 

Fig.  4  shows  the  change  in  the  average  volume  fraction  of  the 
active  material  in  the  spinel  cathode  with  cycle  number  in  the  four 
case  studies.  Initially,  the  volume  fraction  of  the  active  material  in 
the  spinel  cathode  is  0.304.  The  average  volume  fraction  of  the 
active  spinel  cathode  decreases  with  the  cycle  number  due  to  the 
Mn  dissolution  reaction.  Fig.  4  shows  that  the  degradation  rate  of 
the  active  spinel  cathode  is  slow  for  the  cell  cycled  at  the  ambient 
temperature  of  25  °C,  but  the  degradation  of  the  cathode  is 
enhanced  significantly  when  the  ambient  temperature  is  as  high  as 
55  °C.  Fig.  4  also  indicates  that  for  the  same  ambient  temperature, 
the  cell  loses  more  active  material  when  it  is  cycled  within  the 
voltage  range  3.5— 4.3  V  than  the  cell  cycled  within  the  voltage 
range  3.2— 4.0  V  due  to  the  difference  in  the  operating  time.  The 


Fig.  3.  Normalized  discharge  capacity  vs.  cycle  number. 


Fig.  5.  Averaged  film 


(  end  of  discharge  vs.  cycle  number. 
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theoretical  capacity  remaining  in  the  spinel  cathode  in  the  cell 
cycled  with  Tair  —  55  °C  and  voltage  range  3.5— 4.3  V  is  about  76.28%. 

The  average  film  resistance  varying  with  cycle  number  is  pre¬ 
sented  in  Fig.  5  for  the  four  cases.  The  Mn  dissolution  reaction 
which  is  shown  in  Eq.  (1)  not  only  causes  the  loss  of  active  spinel 
cathode  but  also  causes  the  increase  in  the  film  resistance  on  the 
spinel  cathode.  As  the  Mn  dissolution  reaction  proceeds  during 
cycling  of  the  cell,  the  thickness  of  the  inactive  layer,  which  is 
shown  in  Fig.  1  and  assumed  electronically  conductive,  is  growing. 
As  a  result  of  the  growing  inactive  film,  the  film  resistance  increases 
accordingly.  The  temperature  and  voltage  range  effects  on  the  film 
resistance  is  similar  but  in  the  opposite  direction  (increasing  film 
resistance  while  the  volume  fraction  of  the  active  cathode 
decreases)  as  on  the  loss  of  the  active  spinel  cathode  volume 
fraction. 

Fig.  6  shows  the  cell  voltages  during  the  2C  discharge  within 
3.5-4.3  V  and  at  the  ambient  temperature  of  55  °C  in  the  first  cycle 
and  the  last  (50th)  cycle.  The  discharge  capacity  remaining  is  about 
67%  for  the  cell  cycled  under  the  condition:  Tair  =  55  °C  and  voltage 
range  3.5— 4.3  V,  while  the  remaining  theoretical  capacity  in  the  cell 
is  about  76.28%  as  shown  in  Fig.  4.  The  difference  between  the 
available  capacity  and  the  theoretical  capacity  is  due  to  the  differ¬ 
ence  in  the  cell  resistance.  Fig.  6  shows  that  the  potential  drop  from 
the  end  of  charge  (4.3  V)  to  the  beginning  of  a  2C  discharge 
increases  about  1.5  times  in  the  50th  cycle  compared  to  the  first 
cycle. 

Fig.  7  shows  the  normalized  radius  change  with  the  cycle 
number.  The  radii  of  the  active  material  (Ra,p)  and  the  solid  particle 
(Rsp)  are  normalized  by  dividing  by  the  original  radius  of  the 
particle  (Rp)  in  the  spinel  cathode.  The  four  shadows  in  Fig.  7,  each 
of  which  is  enclosed  by  two  thick  lines  (top  line  denotes  the 
normalized  radius  of  the  solid  particle  and  the  bottom  line  denotes 
the  normalized  radius  of  the  active  cathode  part),  are  obtained  from 
the  four  case  studies.  In  each  case,  the  difference  between  the 
radius  of  the  solid  particle  and  the  radius  of  the  active  cathode  part 
indicates  the  thickness  of  the  inactive  film.  The  change  in  the  radius 
of  the  active  part  in  the  spinel  cathode  is  consistent  to  the  change  in 
the  volume  fraction  of  the  active  material  which  is  shown  in  Fig.  4. 
Eqs.  (18)  and  (20)  indicate  that  as  time  goes  infinity,  the  conversion 
of  the  Mn  dissolution  reaction  is  infinity,  as  a  result,  the  volume  of 
the  active  part  is  vanished  and  the  volume  of  the  remaining  inactive 
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Fig.  7.  Radius  change  with  the  cycle  number. 


solid  is  3/4  of  the  initial  volume  of  the  active  cathode.  Consequently, 
as  time  goes  to  infinity  Ra,p  is  zero  and  Rsp  =  v/0.75Rp  which  are 
indicated  in  Eqs.  (21)  and  (22),  respectively. 

The  mechanism  of  Mn  dissolution  mentioned  in  Eqs.  (2)— (4) 
only  considers  the  chemical  reaction  effects.  The  parameter 
values  for  ko  and  Ear  used  in  Eq.  (4)  were  obtained  from  the 
experiments,  [14]  which  consisted  of  immersing  LiMn204  in  the 
electrolyte  at  different  temperatures.  It  has  been  reported  that  the 
Mn  dissolution  rate  also  depends  on  the  state  of  charge,  [6]  the 
composition  of  the  spinel  electrode,  [23]  the  type  of  solute  in  the 
electrolyte,  [24]  etc.  Our  future  work  will  be  to  develop  a  modified 
mechanism  of  Mn  dissolution  by  including  an  electrochemical 
reaction  effect. 


4.  Conclusion 

A  mathematical  model  for  Mn  dissolution  in  a  lithium  ion  cell 
with  a  spinel-based  cathode  was  developed  and  was  applied  in  an 
electrochemical  thermal  model  to  predict  the  cathode  degradation 
during  a  constant  current  charge/discharge  cycling  process  under 
two  selected  ambient  temperatures  of  25  and  55  °C  and  within  two 
selected  voltage  ranges  of  3.2-4.0  V  and  3. 5-4.3  V.  It  was  assumed 
in  the  proposed  model  that  the  cathode  degradation  is  due  to  the 
Mn(IlI)  disproportionation  reaction:  2Mn(III)(s)  -*■  Mn(II)(i)  + 
Mn(IV)(S).  As  a  result,  the  Mn2+  ions  dissolve  in  the  electrolyte 
which  may  transfer  to  the  anode  and  deposit  on  the  anode  to  lead 
to  anode  degradation.  On  the  other  hand,  Mn(IV)  which  is  inactive 
remains  in  the  solid  phase  and  causes  an  increase  in  the  cathode 
resistance.  The  proposed  model  includes  the  effects  of  change  in 
the  radius  of  the  cathode  active  material  particles,  change  in  the 
cathode  porosity  and  change  in  the  film  resistance  on  the  cathode 
degradation.  The  case  studies  reveal  that  the  cathode  loses  its  active 
material  as  the  ambient  temperature  increases  and  the  voltage 
range  goes  higher  where  the  cell  operational  time  is  higher.  The 
case  studies  also  indicate  that  the  ambient  temperature  plays  the 
most  important  role  in  the  Mn  dissolution  reaction. 
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Appendix  A 

The  P2D  model  is  summarized  as  follows. 
In  the  cathode: 


The  solution  phase  concentration  in  the  cathode  is  given  by  the 
material  balance  on  the  electrolyte  phase: 


(Deff4x) + 0  ■  f0)apjp 


The  governing  equation  for  the  charge  balance  in  the  solid  phase 
in  the  cathode: 


3li 


(39) 


where  ap  is  the  specific  area  of  the  active  material  in  the  cathode,  F 
is  Faraday’s  constant,  jp  is  the  pore  wall  flux  and  the  solid  phase 
current  density,  ii,  is  defined  by: 


h 


(40) 


where  <7eff,p  is  the  effective  electronic  conductivity  of  the  cathode, 
and  0i  is  the  solid  phase  potential.  The  effective  electronic 
conductivity  for  the  electrodes  is  determined  by  the  following 
equation: 


where  rp  is  the  porosity  of  the  cathode,  DefflP  is  the  effective  diffu¬ 
sion  coefficient  of  Li+  in  the  binary  electrolyte  in  the  cathode.  The 
initial  concentration  of  the  electrolyte  is  set  to  a  known  constant, 

c0: 


c(t  =  0)  =  c0  (49) 

The  flux  of  Li+  in  the  electrolyte  at  the  cathode/current  collector 
interface  is  zero: 


PX^Lp+Ls+Ln 

The  flux  is  continuous  at  the  separator/cathode  interface: 


(41) 


In  the  separator: 


where  <7;  denotes  the  electronic  conductivity  of  the  dense  material 
in  electrode  i  and  brugg,  is  the  Bruggeman  number  of  electrode  /. 

At  the  cathode/current  collector  interface  (x  =  Lp  +  Ls  +  Ln),  the 
entire  current  density  is  carried  by  the  solid  phase.  That  is: 

k  x=Lp L„  =  lapp  (42) 

where  fapp  is  the  applied  current  density  (the  current  which  is 
divided  by  the  projected  electrode  area),  /app  is  positive  when 
charging  the  cell  and  negative  when  discharging  the  cell. 

At  the  separator/cathode  interface,  no  current  is  carried  by  the 
solid  phase: 

>'lx  L„,u  =  0  (43) 


The  solution  phase  concentration  in  the  separator  is  governed 
by: 


where  cs  is  the  porosity  of  the  separator  and  Det f,s  is  the  diffusion 
coefficient  of  Li+  in  the  electrolyte  in  the  separator.  The  initial 
concentration  of  the  electrolyte  is  set  to  a  known  constant,  Co: 

c(t  =  0)  =  c0  (53) 

The  fluxes  of  Li+  in  the  electrolyte  at  the  anode/separator  and 
separator/cathode  interfaces  are  continuous: 


The  charge  balance  in  the  electrolyte  is  given  by: 


di2 

dx 


OpFjp 


(44) 


where  the  current  density  in  the  electrolyte  filled  in  the  cathode,  t'2, 
is  defined  by: 

c«) 

where  xeff,p  is  the  effective  ionic  conductivity  of  the  electrolyte  in 
the  cathode,  0 2  is  the  electrolyte  potential,  R  is  the  universal  gas 
constant,  T  is  the  temperature  and  tp  is  the  transference  number 
and  c  is  the  concentration  of  the  binary  electrolyte. 

At  x  =  Lp  +  Ls  +  Ln,  there  is  no  current  in  the  electrolyte  passing 
through  the  cathode/current  collect  interface: 

‘2x  -Lp  -Is ;  l.„  =  0  (46) 

At  x  =  Ln  +  Lj,  the  current  density  in  the  electrolyte  at  the 
separator/cathode  interface  is  continuous: 

i2lx=(L„+Ls)-  =  >2lx=(Ln+Ls)+  (47) 


— °eff.n  S 


_  =  —  Deff.s7 


(54) 


-Deff,s 


0C| 

S0x|, 


x=(Ln+Lsr 


del 

p9x|j 


■£=(L„+Ls)+ 


(55) 


where  Deff,n  is  the  diffusion  coefficient  of  Li+  in  the  electrolyte  in  the 
anode.  It  should  note  that  the  concentration  of  Li+  at  the  interfaces 
is  also  continuous,  in  another  words,  there  is  no  concentration 
jump  at  the  interfaces. 

Because  there  is  no  active  solid  phase  in  the  separator,  all 
currents  are  carried  by  the  electrolyte,  that  is: 


where  the  current  density  in  the  electrolyte  in  the  separator  is 
given  by: 


,  „  602  2xeffsRT  /  o  \  dine 


(57) 


where  Keff,s  is  the  effective  ionic  conductivity  of  the  electrolyte  in 
the  separator. 
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The  current  densities  at  the  anode/separator  and  separator/ 
cathode  interfaces  are  continuous: 


where  Ds  n  is  the  diffusion  coefficient  of  Li+  in  the  particles  in  the 
anode.  The  initial  condition  is: 


!2|*=l„  =  h\x=K  (58) 

MaHM-L*)*  =  ,2lx=(Ln+u)+  (59) 


Cs.n(t  =  0)  =  Cs  n  o  (70) 

where  cs>n,o  is  the  initial  concentration  of  Li+  in  the  particles  in  the 
anode.  And  the  boundary  conditions  for  the  anode  particles  are: 


In  the  anode: 


The  charge  balance  in  the  solid  phase  in  the  anode  is  given  by: 


where  an  is  the  specific  of  the  active  material  in  the  anode,  j„  is  the 
pore  wall  flux  in  the  anode  and  the  current  density  in  the  solid 
phase  in  the  anode,  i'i,  is  defined  by: 


where  Rn  is  the  radius  of  the  particles  in  the  anode. 


(71) 


(72) 


it 


(61) 


where  <7eff,n  is  the  electronic  conductivity  of  the  solid  phase  in  the 
anode.  The  boundary  conditions  for  fa  are: 

iilx=o  =  lapp  and  ii  |x=Ln  =  0  (62) 


Appendix  B 

Expressions 

The  pore  wall  flux  at  the  surface  of  the  particles  in  electrode  i  is 
related  to  the  potentials  in  the  solid  and  solution  phases  through 
the  Butler— Volmer  type  kinetics  equation: 


The  potential  in  the  solution  phase  in  the  anode  is  given  by: 

(63) 

where  the  current  density  in  the  electrolyte  in  the  anode,  h,  is 
defined  by: 


902  ,  2iceS  nRT  f 


6lnc 


(64) 


where  /ceff,n  is  the  effective  ionic  conductivity  of  the  electrolyte  in 
the  anode. 

The  boundary  conditions  for  (j>2  in  the  anode  are: 

02lx=o  =  0  and  i2|x=Ln  =  t2|x=I+  (65) 


where  the  solution  phase  potential  at  x  =  0  is  set  as  the  reference 
potential. 

The  material  balance  in  the  electrolyte  in  the  anode  is  given  by: 


9(fnc) 

at 


(D«»S)  +  ('-SW"  1661 


where  rn  is  the  porosity  of  the  anode,  Deff,n  is  the  diffusion  coeffi¬ 
cient  of  Li+  in  the  binary  electrolyte  in  the  anode.  The  initial 
concentration  of  the  electrolyte  is  set  to  a  known  constant,  cp: 


where  kj  is  the  reaction  rate  constant  in  electrode  i,  cs,;iSurf  is  the 
surface  concentration  of  Li+  in  the  particles  in  electrode  i,  cs,i,max  is 
the  maximum  concentration  of  Li+  in  the  particles  in  electrode  i,  rji 
is  the  overpotential  on  electrode  i.  The  overpotentials  on  the 
cathode  and  the  anode  are  defined  as  follows: 

Vp  =  0i, p  -  02, p  ~  Up  ~  RfnmFjp  (74) 

%  =  01, n  -  02, n  -  Un  (75) 

The  open  circuit  potentials  in  the  electrodes  are  functions  of 
temperature  and  states  of  charge  (SOC)  of  the  electrodes  as  follows: 

Ui  =  Uj  ref  +  (T|  —  Tref)  (76) 

where  l/;,ref  is  the  open  circuit  potential  as  a  function  of  state  of 
charge  (SOC)  measured  at  the  reference  temperature,  Tref. 

The  open  circuit  potentials  at  Tle f  =  25  °C  for  the  LiMn204 
cathode  and  for  the  carbon  anode  as  functions  of  state  of  charge  are 
given  by  [20]: 


c(t  =  0)  =  c0  (67) 

The  flux  of  Li+  at  the  current  collector/anode  interface  is  zero: 

-d«4L = 0  <68) 

The  flux  at  the  anode/separator  interface  is  continuous  as  shown 
in  Eq.  (54). 

The  concentration  of  Li+  in  the  particles  in  the  anode  is  given  by 
the  Fick’s  law  in  the  spherical  coordinates: 


Up, ref  =4.19829  +  0.0565661  tan  h(-14.55460p  +  8.60942) 
-  0.0275479  ( 


\(0. 998432  -  i9p)u-4^4“ 
-  O.157123exp(-O.O47380p) 

4-  0.810239exp(  -  4O(0P  -  0.133875)) 


and 


9cs,n  =  1  9/ 2D  9Cs,n\  f691  Un  ref  =  -0.16  +  1.32exp(-30n)  +  lOexp(-2OOO0n) 

5 1  r2dr{  s,n  9r  J  K  ’ 


(78) 
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where  the  SOC  is  defined  by: 

0i  =  §7™x’  l  =  p,n  (79) 

The  [aU/9TL  data  as  a  function  of  SOC  for  LiMn204  reported  by 
Thomas  et  al.  [25]  was  refitted  to  an  empirical  equation  presented 
by  Srinivasan  and  Wang  [26]  as  follows: 


where  DSi,-|[ef  is  the  diffusion  coefficient  of  Ii+  in  the  solid  particle  in 
electrode  i,  Efld,i  is  the  activation  energy  related  to  the  diffusion 
coefficient  of  Li+  in  the  solid  phase  in  electrode  i. 

The  specific  area  of  the  active  material  in  electrode  i  is  related  to 
the  volume  fraction  of  the  active  materials  in  electrode  i  as  follows: 

aP  =  ^  +p  and  a„  =  ^-c„,n  (84) 


[w]  =  a'exP(a2^p)  +  <23Sin(a40p)  +  a5sin(a6tfp  +  a7) 

+  a8sin(ag0p  +  am)  +  bo  +  b\6p  +  b20p  +  b3$p 
+  b4(exp(^-9p^b5yj  (80) 


The  concentration  and  temperature  dependent  ionic  conduc¬ 
tivity  and  the  diffusion  coefficient  in  the  binary  electrolyte  are 
given  by  [27]: 


(  -10.5  +  0.668  x  10  3c,-  +  0.494  x  10“6c?  +  0.0747) 
y  -1.78  x  10-5CjTi  -  8.86  x  1010c?7)  -  6.96  x  lO  5!?  +  2.80  x  10  8c,r? 


(85) 


where:  ai  =  -272.45547143156853,  a2  =  -8.4660463364774969, 
a3  =  3.9548818181675394,  a4  =  -6.3644778205164450, 

a5  =  -0.11494477221293507,  a6  =  -26.341590055374475, 
a7  =  28.202482124134448,  a8  =  -0.0045805692645594738, 
a9  =  114.74983936455860,  aw  =  -72.387775508734848, 

b0  =  53.722802544438963,  fa,  =  -164.19140744090566, 
fa2  =  136.83084104679463,  b3  =  -26.071652680281360, 

fa4  =  0.74302471734078734,  b5  =  -0.53278657307359123, 
b6  =  0.18367552820218680,  and  the  unit  for  [8U/3T]p  is  mV  K"1. 

The  [dU/dT]n  as  a  function  of  SOC  for  a  carbon-based  anode 
(meso  carbon  micro  beads,  MCMB)  was  measured  at  USC  and  was 
fitted  to  the  following  empirical  equation: 

rain  _  ao  +  a^n  +  ^n  +  ^n  +  ^n  +  ^n  +  ^n  +  ^n  +  ^n 
1+bA  +  Mn  +  Mn  +  Mn  +  Mn  +  Mn  +  Mn  +  Mn 
(81) 

where:  a0  =  0.005269056,  ai  =  3.299265709000005, 

a2  =  -91.79325798000001,  a3  =  1004.911008,  a4  =  -5812.278127, 
a5  =  19,329.7549,  a6  =  -37,147.8947,  a7  =  38,379.18127, 
as  =  -16,515.05308,  b i  =  -48.09287227,  b2  =  1017.234804, 
b3  =  -10,481.80419,  b4  =  59,431.30001,  b5  =  -195,881.6488, 
be  =  374,577.3152,  b7  =  -385,821.1607,  b8  =  165,705.8597,  and  the 
unit  for  [ai//8T]n  is  mV  K_1. 

The  temperature  dependent  reaction  rate  constants  for  the 
electrodes  are  determined  by: 

fc,  _  t,,«rexp(%£  (,+-  - 1) )  (82) 

where  k,iref  is  the  reaction  rate  constant  for  electrode  i  at  the 
reference  temperature,  Tref,  Eak,i  is  the  activation  energy  related  to 
the  reaction  rate  constant  for  electrode  i. 

The  temperature  dependent  diffusion  coefficient  of  Li+  in  the 
solid  phase  is  given  by: 

(83) 


D,  =  10  4  x  10  4  43-'(,  -  S«i5  3cr°.22x10  3Cj  (86) 

The  effective  ionic  conductivity  and  diffusion  coefficient  in  the 
binary  electrolyte  are  determined  by  the  following  equations, 
respectively: 

Keff,i  =  KiE|>rUSSi!  i  =  Pi  S,  n  (87) 

Deffj  =  D^mm,  i  =  p,  s,  n  (88) 

The  cell  voltage  is  determined  as  the  difference  in  the  solid 
phase  potentials  between  the  two  ends  of  the  cell: 

^cell  =  01  lx=L„+I,+Lp  -  0llx=O  (89) 


List  of  symbols 


a,-  specific  surface  area  of  electrode  i  (i  =  p,  n),  m2  m  3 

bruggj  Bruggeman  coefficient  of  region  i  (i  =  p,  s,  n) 
c  electrolyte  concentration,  mol  m-3 

Co  initial  electrolyte  concentration,  mol  m-3 

cs,i  concentration  of  lithium  ions  in  the  intercalation  particle 
of  electrode  i  ( i  =  p,  n),  mol  m  3 
cSijro  initial  concentration  of  lithium  ions  in  the  intercalation 
particle  of  electrode  i  (i  =  p,  n),  mol  m  3 
cs,i,max  maximum  concentration  of  lithium  ions  in  the 

intercalation  particle  of  electrode  i  (i  =  p,  n),  mol  m  3 
D  electrolyte  diffusion  coefficient,  m2  s_1 

DSj  lithium  ion  diffusion  coefficient  in  the  intercalation 
particle  of  electrode  i  (i  =  p,  n),  m2  s_1 
F  Faraday’s  constant,  96,487  C  equiv  1 

bcoef  heat  transfer  coefficient,  W  rrT2  K-1 
Japp  applied  current  density,  A  m  2 

11  solid  phase  current  density,  A  m-2 

12  solution  phase  current  density,  A  m”2 
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wall  flux  of  Li+  on  the  intercalation  particle  of  electrode  i 
(t  =  p,  n),  mol  m-2  s_1 

intercalation/deintercalation  reaction  rate  constant  of 

electrode  i  ( i  =  p,  n),  mol0  5  m2-5  s-1 

thickness  of  region  i  (i  =  p,  s,  n),  m 

negative  electrode 

positive  electrode 

radial  coordinate,  m 

universal  gas  constant,  J  mol-1  K-1 

radius  of  the  active  spinel,  m 

film  resistance  at  the  cathode,  O  m2 

initial  radius  of  the  intercalation  particle  of  electrode  i 

(i  =  p,  n),  m 

radius  of  the  solid  particle  in  the  cathode,  m 
Li+  transference  number  in  the  electrolyte 
absolute  temperature,  K 

open  circuit  potential  of  electrode  i  ( i  =  p,  n),  V 
spatial  coordinate,  m 
porosity  of  region  i  ( i  =  p,  s,  n) 
volume  fraction  of  fillers  of  electrode  i  (i  =  p,  n) 
volume  fraction  of  the  inactive  material  in  cathode 
dimensionless  concentration  of  lithium  ions  in  the 
intercalation  particle  of  electrode  i  (ft  =  cs,i, surf/cs,i,max) 
ionic  conductivity  of  the  electrolyte,  S  m-1 
effective  ionic  conductivity  of  the  electrolyte  in  region  i 
(t  =  p,  s,  n),  S  m'1 

electronic  conductivity  of  the  solid  phase  of  electrode  i 
(i  =  p,  n),  S  m  1 

effective  electronic  conductivity  of  the  solid  phase  of 
electrode  i  (i  =  p,  n),  S  m-1 
solid  phase  potential,  V 
solution  phase  potential,  V 


References 


[1]  M.  Wohlfahrt-Mehrens,  C.  Vogler,  J.  Garche,  J.  Power  Sources  127  (2004) 
58-64. 

[2]  W.  Choi,  A.  Manthiram,  J.  Electrochem.  Soc.  153  (2006)  A1760-A1764. 

[3]  P.  Arora,  R.E.  White,  M.  Doyle,  J.  Electrochem.  Soc.  145  (1998)  3647-3667. 

[4]  Y.Y.  Xia,  Y.H.  Zhou,  M.  Yoshio,  J.  Electrochem.  Soc.  144  (1997)  2593-2600. 

[5]  D.H.  Jang,  Y.J.  Shin,  S.M.  Oh,  J.  Electrochem.  Soc.  143  (1996)  2204-2211. 

[6]  L.F.  Wang,  C.C.  Ou,  K.A.  Striebel,  J.J.S.  Chen,  J.  Electrochem.  Soc.  150  (2003) 
A905-A911. 

[7]  T.  Aoshima,  K.  Okahara,  C.  Kiyohara,  K.  Shizuka,  J.  Power  Sources  97  (2001) 
377-380. 

[8]  T.  Inoue,  M.  Sano,  J.  Electrochem.  Soc.  145  (1998)  3704-3707. 

[9]  D.  Aurbach,  M.D.  Levi,  K.  Gamulski,  B.  Markovsky,  G.  Salitra,  E.  Levi,  U.  Heider, 
L.  Heider,  R.  Oesten,  J.  Power  Sources  81  (1999)  472-479. 

[10]  D.  Kim,  S.  Park,  O.B.  Chae,  J.H.  Ryu,  Y.U.  Kim,  R.Z.  Yin,  S.M.  Oh,  J.  Electrochem. 
Soc.  159  (2012)  A193-A197. 

[11]  T.  Eriksson,  A.M.  Andersson,  A.G.  Bishop,  C.  Gejke,  T.  Gustafsson,  J.O.  Thomas, 
J.  Electrochem.  Soc.  149  (2002)  A69-A78. 

[12]  D.  Aurbach,  K.  Gamolsky,  B.  Markovsky,  G.  Salitra,  Y.  Gofer,  U.  Heider, 
R,  Oesten,  M.  Schmidt,  J.  Electrochem.  Soc.  147  (2000)  1322-1331. 

[13]  J.  Park,  J.H.  Seo,  G.  Plett,  W.  Lu,  A.M.  Sastry,  Electrochem.  Solid-State  Lett.  14 
(2011)  A14-A18. 

[14]  C.H.  Lu,  S.W.  Lin,  J.  Mater.  Res.  17  (2002)  1476-1481. 

[15]  M.  Doyle,  T.F.  Fuller,  J.  Newman,  J.  Electrochem.  Soc.  140  (1993)  1526-1533. 

[16]  V.  Srinivasan,  J.  Newman,  J.  Electrochem.  Soc.  151  (2004)  A1517-A1529. 

[17]  Q.  Zhang,  R.E.  White,  J.  Electrochem.  Soc.  154  (2007)  A587-A596. 

[18]  B.C.  Han,  A.  Van  der  Ven,  D.  Morgan,  G.  Ceder,  Electrochim.  Acta  49  (2004) 
4691-4699. 

[19]  L.  Cai,  R.E.  White,  J.  Electrochem.  Soc.  157  (2010)  A1188-A1195. 

[20]  M.  Doyle,  J.  Newman,  A.S.  Gozdz,  C.N.  Schmutz,  J.M.  Tarascon,  J.  Electrochem. 
Soc.  143  (1996)  1890-1903. 

[21]  SUNDIALS,  online:  https://computation.llnl.gov/casc/sundials/main.html, 
(accessed  20.03.12). 

[22]  COMSOL  Multiphysics®,  online:  http://www.comsol.com/,  (accessed  20.03.12). 

[23]  R.J.  Gummow,  A.  de  Kock,  M.M.  Thacheray,  Solid  State  Ionics  69  (1994)  59-67. 

[24]  A.  Du  Pasquier,  A.  Blyr,  P.  Courjal,  D.  Larcher,  G.  Amatucci,  B.  Gerand, 
J.-M.  Tarascon,  J.  Electrochem.  Soc.  146  (1999)  428-436. 

[25]  ICE.  Thomas,  C.  Bogatu, J.  Newman, J.  Electrochem.  Soc.  148  (2001 )  A570-A575. 

[26]  V.  Srinivasan,  C.Y.  Wang,  J.  Electrochem.  Soc.  150  (2003)  A98-A106. 

[27]  L.O.  Valoen,  J.N.  Reimers.J.  Electrochem.  Soc.  152  (2005)  A882-A891. 


